Optimal laser-control of double quantum dots 
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Coherent single-electron control in a realistic semiconductor double quantum dot is studied the- 
oretically. Using optimal-control theory we show that the energy spectrum of a two-dimensional 
double quantum dot has a fully controllable transition line. We find that optimized picosecond 
laser pulses generate population transfer at significantly higher fidelities (> 0.99) than conventional 
sinusoidal pulses. Finally we design a robust and fast charge switch driven by optimal pulses that 
are within reach of terahertz laser technology. 



Double quantum dots (DQDs), i.e., coupled two- 
dimensional (2D) electron traps, have been under recent 
and extensive studies both experimentally^^^^^ and theo- 
retically.^^^^^ The main interest in DQDs arises from their 
potential for solid-state quantum computation that could 
be achieved in principle by rapidly switching voltages of 
electrostatic gates. The gates permit to tune at will the 
system geometry and hence the electronic properties of 
DQDs. Coherent manipulation of a single charge^ and 
coupled spins^ has already been achieved, and recently a 
coherence time of ~ 200 ns was obtained for a well iso- 
lated silicon DQD.^ Theoretical studies on single-electron 
transport inside the DQD driven by linear switches and 
linearly polarized continuous waves (CWs) were reported 
very recently.^ In the latter case the transport is rather 
sensitive to possible anharmonicity of the potential and 
limited to uncoupled dots far apart from each other. 
Electron control in DQDs has been studied also using 
genetic algorithms^ as well as rotating-wave and reso- 
nant approximations leading to a reduction to a three- 
level system^^ To the best of our knowledge, however, a 
general A^-level control scheme by using direct external 
electric fields has not been introduced for 2D-DQDs until 
now. 

In this paper we discuss the controllability criteria for 
single-electron states of DQDs by means of external laser 
pulses. We show that at certain interdot distances some 
of the single-electron states allow full population transfer 
from the ground state to those states. We apply quantum 
optimal-control theory (OCT)^ which yields the optimal 
laser pulses for predefined transitions. We obtain high 
occupations (> 99%) of the target states in a realistic 
DQD in a few picoseconds, which is well in the coher- 
ent regime. If the initial and final states are chosen to 
have full localization of the electron in one or the other 
dot, this scheme enables rapid and controlled transport 
which is not sensitive to the interdot distance or to the 
inevitable anharmonicities in the confining potential. 

In the static 2D Hamiltonian, Hq = — [d^ + 9^) /2 + 
Vc{x^ y), the external potential describing the DQD is, in 
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FIG. 1: (color online). Left panel: Lowest eigenenergies of 
a double quantum dot with ooq = 0.5 as a function of the 
interdot distance. Black, red (thick), and blue (thin) curves 
mark the ground state, controllable states, and uncontrollable 
states, respectively. Right panel: Densities of six lowest eigen- 
states at d = 5. The dashed lines mark the nodes of the wave 
functions. 



its most common form,^ given by 



Vc{x,y) ^min 



, (1) 



where d is the distance between the potential minima, 
and coq = 0.5 is the confinement strength with a typical 
value for DQDs. We apply the effective-mass approxima- 
tion for electrons moving in GaAs and set the effective 
mass to m* = O.OOTmg and the dielectric constant to 
hi = 12.7eo. The energies, lengths, and times are given in 
effective atomic units (a.u.): Ha* = (m*/A>:^)Ha ^ 11.30 
meV, ttg = {m^ / f<i)ao ~ 10.03 nm, and = fi/Ha* ~ 
58.23 fs, respectively. 

In Fig. [1] we show the lowest energy levels as a func- 
tion of d (left panel) as well as densities of six lowest 
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eigenstates at d = 5 (right panel). The d = limit 
corresponds to the well-known shell structure of a single 
2D harmonic oscillator (HO), where the energy levels are 
n-fold degenerate at energies nuoo {n = 1,2,...). Increas- 
ing d leads to lifting of the degeneracies and generates 
crossings and avoided crossings between the energy lev- 
els. In the weak-coupling limit d ^ oo the bundling of 
the levels at ncoo is restored. In this limit, the energies 
are 2n-fold degenerate corresponding to two (uncoupled) 
HOs. We label the states in the DQD as IV^ij) = 
where i = 0, 1, . . . denotes the {i + l)th bundle of states 
at (i = 0, and j = 0, 1, . . . denotes the {j + l)th state in 
each bundle. As visualized in the right panel of Fig. [U 
index i also corresponds to the number of nodes in the 
wave function. We show below that our labeling yields 
simple rules for controllable and uncontrollable states as 
a function of d. 

It is a well-known fact in control theory that an infinite- 
level single HO is not controllable in the dipole approx- 
imation^i^ This is in contrast with the truncated HO 
which (in most cases) satisfies the controllability crite- 
ria^ii In qualitative arguments, the uncontrollability of 
a single HO {d = 0) in the dipole approximation stems 
from the equidistant single-electron level spacings. Using 
the above labeling of the states, the nonvanishing dipole- 
matrix elements {ij\jl\kl) between the HO states {d = 0) 
are 
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where fi = — er is the dipole operator. Eq. ([2]) holds 
for all values of d, and the corresponding energy-level 
spacings remain constant. Hence, transitions ij ^ (i ± 
l)(j ± 1) remain uncontrollable. On the other hand, the 
elements given in the LHS of Eq. (j3|) change as a func- 
tion of as well as the level spacings. The behavior 
already suggests that transitions ij ^ (i ± l)j become 
controllable when d is increased from zero. Our calcu- 
lations below confirm this prediction. Starting with the 
ground state |00), the controllable transition line is then 
|00) |10) |20) . . . (see the red curves in Fig. [1]). In 
the weak-coupling limit (d oo), however, the LHS of 
Eq. (j3|) becomes 



lim I 

d^oo 



OO if i — j is even 

V^i^- j) if ^ - J is odd. 



(4) 

so that only the trivial transitions between the degen- 
erate gerade and ungerade states are possible. In the 
large-d regime, however, breaking the inter dot symmetry 
leads to the possibility of charge transport between the 
dots (see below). 

Next we apply OCT in order to find optimal laser 
pulses for transitions from the initial state |$i) = |^(t = 
0)) to the target state |<I>f) = |^(t = T)) in a fixed 
time interval T. In the OCT formalism we maximize the 
overlap | (^(T)|<I>f) P while minimizing the fiuence of the 
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FIG. 2: (color online). Maximum occupation (logarithmic 
scale) of the target state as a function of the interdot distance 
d in transitions |00) |10) with pulse length T = 50 (circles), 
and 1 00) |20) with T = 100 (squares). The lines are to 
guide the eye. The dashed line denotes the maximum target- 
state occupation 0.6 for a single harmonic oscillator. The 
jump marked by an arrow is due to a resonance effect (see 
text). 



laser pulse. The control equations ar( 



.12 



e(t) 



H^{t), ^(0) = $1, 



a 



-Im {x{t)m{t)), 



(5) 
(6) 

(7) 



where Eq. ([5]) is the time-dependent Schrodinger equa- 
tion with H = Hq — /i6(t), and x{t) is the Lagrange 
multiplier. The optimal laser pulse e{t) is provided at 
the end of the iterative procedure. We point out that 
the initial pulse (zeroth iteration) is sinusoidal and has 
both X and y components, whereas the converged opti- 
mal pulse is always found to be polarized in the x di- 
rection, i.e., e{t) = e{t)x. In Eq. ([7j) we choose a si- 
nusoidal pulse envelope A{t) = sin^(7rt/T), and restrict 
the pulse intensity by a penalty factor a. Unless stated 
otherwise, we have fixed a = 0.5 leading to pulse inten- 
sities 10^ ... 10^ W/cm^. We apply a rapidly converging 
numerical schemei2.Lll that has been implemented in the 
octopus code^i^ 

We consider excitations only from the ground state and 
set l^i) = I go). Figure [2] shows the maximum overlaps 
I (^(T)I^f) i-e., maximum occupations of the target 
states |10) and |20) as a function of the interdot distance 
(note the logarithmic scale). The pulse lengths are fixed 
to T = 50 and T = 100, respectively. When d is in- 
creased, the target-state occupations increase from the 
HO value < 0.6 (see Ref. [l^) exponentially to close to 
one. As expected, at large interdot distances {d > 8), 
corresponding to the uncoupling of the DQD, the occu- 
pations for |00) |20) decrease back to the HO value 
marked by a dashed line. On the other hand, the oc- 
cupation for 1 00) ^ 1 10) increases even further in this 
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FIG. 3: (color online). Upper panel: Optimized pulses (x 
components) for transitions |00) |10) (a) and |00) |20) 
(b). The interdot distances are fixed to d = 3 and 5 and the 
pulse lengths to T = 50 and 100, respectively. Lower panel: 
Occupations of states involved in the transitions. 

limit due to the asymptotic degeneracy of the states (see 
above). Generally, for this transition a pulse length of 
T = 50 a.u. ~ 3 ps is sufficient to achieve high occupa- 
tions. In the transition |00) |20) instead, T = 100 a.u. 
~ 6 ps is needed for the same accuracy. 

As seen in Fig. [21 the target-state occupations above 
90% are generally obtained at distances 2 < d < 8, i.e., 
at (i ~ 20 ... 80 nm. This length scale is well realizable in 
experiments. ^'^'^ We point out that the most distinctive 
deviations in the occupations in this regime arise from 
resonance effects. For example, the jump in the |00) ^ 
|10) transition at (i = 3.6 . . .3.7 (see the arrow in Fig. [2]), 
where the occupation decreases from 0.9997 to 0.9591, 
is due to the degeneracy of cJqq and ujIq — cjfo- The 
degeneracy disturbs the optimal transition path leading 
to reduced maximum occupation. 

In Fig. [3] we give two examples of optimized pulses and 
the occupations of states during transitions |00) |10) 
when d = 3 (a,c) and |00) ^ |20) when d = b (b,d). The 
pulse for |00) |10) has the resonant frequency cJqo 
the major component, but as seen in the occupations, 
state |20) is also populated during the transition. For 
comparison, we exposed the system also to a CW having 
the resonant (Rabi) frequency Qr = cJqq and the same 
fixed length (T = 50) as the optimized pulse. We set the 
pulse amplitude to flR/fiQQ = t:/{ii^T) to satisfy the 
TT-pulse condition. This leads to occupation 0.93 of state 
|10>. In view of the 0.999 occupation achieved by the 
optimal pulse, OCT is superior to the CW approach. A 
comprehensive comparison of occupations achieved with 
different pulse lengths applying OCT and CWs for a 
quantum-ring system is given in Ref. [12]. 

As shown in Fig.[3Kd), the |00) — > |20) transition must 
be mediated by the |10) state since a direct transition 



is forbidden by the dipole selection rules. However, the 
intermediate state |10) does not need to get fully popu- 
lated during the optimized transition process. The higher 
states are also involved in the process, and in this exam- 
ple the state |30) receives considerable occupation with a 
maximum of ~ 20% at t ~ 90. The final target-state oc- 
cupation is 0.998, whereas using two resonant CWs one 
after the other we could not exceed 0.9 (pulse length fixed 
to T = 100). Generally, in multilevel transitions of this 
type, the efficiency of OCT is pronounced with respect 
to CWs due to the multiplication of errors in the latter 
approach when full population of intermediate states is 
required before changing the resonant frequency. 

The control of electron transport in the DQD requires 
initialization of the state by localizing the electron in the 
one of the dots. Namely, the single-electron ground state 
of the DQD is a gerade state with half of the electron in 
one well and half in the other (see the lowermost figure 
in the right panel of Fig. [1]). Only in the limit d ^ oo 
the localized states become degenerate eigenstates. How- 
ever, the initial state can be fixed with certainty by ap- 
plying a constant external field or by adjusting the gates 
in the DQD device in order to create a potential shift 
between the two dots. After the initial localization, elec- 
tron transport into another well can be driven using a 
linear switch or a CW with the resonant confinement 
frequency uoq. These types of transport on time scales 
of 10~^^ s, close to the required times for SWAP oper- 
ations in experiments,^ were reported very recently by 
F0rre and co-workers.^ 

Now we show that OCT provides a very fast (switch- 
ing times of a few picoseconds) and stable alternative to 
control the electron transport in a DQD. First we break 
the degeneracy of the ground state by setting the exter- 
nal potential in the (lower) left dot yj^^* - 0.2. 
Then the ground state |00l) = \L) and the first excited 
state |00i^) = \R) correspond to electron localization in 
the (lower) left and (upper) right dot, respectively, pro- 
vided that d is sufficiently large. The result of the OCT 
calculation for transition |L) is shown in Fig. [H 
The pulse length is fixed to T = 100 (~ 6 ps), and the 
interdot distance is d = 6 corresponding to relatively 
large coupling between the dots. The spectrum of the 
optimized pulse (a) has a large peak at the resonant fre- 
quency 00 = 00^ = 0.5 a.u. ~ 9 THz, and a few smaller 
peaks at cc; ~ 0.2... 0.3 and uo 1. The small peaks 
correspond to transitions in the higher states which get 
significantly populated in the transitions. This is visu- 
alized in Fig. [31(b) showing the occupations of |L) and 
\R) (solid lines) that sum up only to about 50% in the 
middle of the transition at T ~ 50. We also plot the 
integrated electron densities pR and pL (dashed lines) 
in the (lower) left and (upper) right parts of the DQD, 
respectively. The quasi-periodic oscillations in the den- 
sities indicate that the electron charge transfers in blobs 
as visualized in the snapshots in Figs.HJc-h). In this ex- 
ample we find the final occupation | {^(T)\R) p = 0.985. 
It is worth noting that the final occupations are closer 
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FIG. 4: (a) Spectrum of the optimized pulse (inset) for the 
electron transport process |L) \R) in a fixed time T = 100. 
The pulse has a rectangular envelope A{t) = 1 and the penalty 
factor is a = 1. (b) Occupations of states \R) and |L) (solid 
lines) and the integrated densities pn and pL (dashed lines) in 
the right and left dot during the transport, (c-h) Snapshots of 
the total electron density pr + pl- The double quantum dot 
has d = 6, a;o = 0.5, and well-depth asymmetry of Vb = 0.2. 

to one at larger interdot distances so that the coupling 
between the dots is weaker. In this regime the optimized 
pulses may attain a linear slope. This immediately sug- 
gests that a combination of a linear field and OCT could 
be ideal in controlling transport in the uncoupled regime. 
In this work, however, we focus on coupled DQDs which 
we in fact find the most challenging in terms of electron 
control in nanoscale applications. 

There are three significant advantages in OCT with re- 
spect to a linear switch or a resonant CW when generat- 
ing electron transport in DQDs. First, the optimization 
procedure is insensitive to the interdot coupling. The 
CW approach in principle requires an uncoupled system 
for electron transport with unit probability."^ We tested a 
CW with the resonant coq frequency for the above exam- 
ple with various amplitudes and could not exceed occupa- 
tion 0.4 of the target state \R). Secondly, the optimized 
pulses can be very short in order to achieve sufficient oc- 



cupation, since the higher-lying states are incorporated 
in the control problem. In the above example the pulse 
length of ~ 6 ps is more than by a factor of ten smaller 
than in linear and CW switches and well below the re- 
cently measured coherence times of hundreds of nanosec- 
onds in DQDs.^ Thirdly, and probably most importantly, 
the OCT approach is insensitive to deviations in the ex- 
ternal potential, since the specific shape of the DQD is 
taken into account explicitly through the external po- 
tential in the Hamiltonian. We tested this by adding a 
fourth-order anharmonic term in the external potential 
[Eq. ([T|)]. This results in a change in the optimal pulse 
shape, and in an increase in the required pulse length 
(in the case of a positive anharmonicity ) . However, we 
found no decrease in the obtained target-state occupa- 
tion when using anharmonicities that lead to dramatic 
loss of accuracy in the CW approach.^ In practical ap- 
plications the exact shape of the external potential could 
be obtained by measuring the differential conductance 
in a single-electron transport experiment, and thereafter 
numerically solving the (inverse) Schrodinger equation. 
This approach has been applied to modeling external im- 
purities inside single quantum dotsJ^ 

We expect that the experimental creation of opti- 
mal laser pulses such as presented in this work will be 
soon within reach of laser technologies. The intensities 
10^ . . . 10^ W/cm^ required for DQD excitations can be 
already obtained in the THz regime by high-power free- 
electron lasers. On the other hand, shaping of picosec- 
ond THz pulses has been recently made possible by em- 
ploying transient polarization grating.^- Quantum cas- 
cade lasers may also provide an applicable route for 
precise pulse shaping in the THz regime. The rapid de- 
velopments in the THz laser technology^- will eventually 
lead to the combination of sufficient pulse power and ac- 
curate manipulation of the pulse shape. 

Finally we point out that in our future work we aim at 
combining our approach with magnetic-field optimization 
which could allow us to coherently control the spin state 
simultaneously with the electron localization. 
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